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Abstract 

During the last years, through the combined effort of the insight, coming from physical 
intuition and computer simulation, and the exploitation of rigorous mathematical methods, 
the main features of the mean field Sherrington-Kirkpatrick spin glass model have been 
firmly established. In particular, it has been possible to prove the existence and uniqueness 
of the infinite volume limit for the free energy, and its Parisi expression, in terms of a 
variational principle, involving a functional order parameter. Even the expected property of 
ultrametricity, for the infinite volume states, seems to be near to a complete proof. 
The main structural feature of this model, and related models, is the deep phenomenon 
of spontaneous replica symmetry breaking (RSB), discovered by Parisi many years ago. By 
expanding on our previous work, the aim of this paper is to investigate a general frame, where 
replica symmetry breaking is embedded in a kind of mechanical scheme of the Hamilton- 
Jacobi type. Here, the analog of the "time" variable is a parameter characterizing the strength 
of the interaction, while the "space" variables rule out quantitatively the broken replica 
symmetry pattern. Starting from the simple cases, where annealing is assumed, or replica 
symmetry, we build up a progression of dynamical systems, with an increasing number of 
space variables, which allow to weaken the effect of the potential in the Hamilton-Jacobi 
equation, as the level of symmetry braking is increased. 

This new machinery allows to work out mechanically the general if-step RSB solutions, in a 
different interpretation with respect to the replica trick, and lightens easily their properties 
as existence or uniqueness. 

1 Introduction 

In the past twenty years the statistical mechanics of disordered systems earned an always in- 
creasing weight as a powerful framework by which analyze the world of complexity [5] [6] 

HU pH] [29] |3I]. 
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The basic model of this field of research is the Sherrington-Kirkpatrick model [26j (SK) for a spin 
glass, on which several method of investigation have been tested along these years [3 J [4J [9j [13J 
|17| p3j [24j |36| |37| . The first method developed has been the replica trick [27j[32j which, in a 
nutshell, consists in expressing the quenched average of the logarithm of the partition function 
Z(/3) in the form E In Z(/3) = lim n ^o E(Z(/3) n — 1) jn. Since the averages are easily calculated for 
integer values of n, the problem is to find the right analytic continuation allowing, in some way, 
to evaluate the n — )■ limit, at the least for the case of large systems. Such analytic continuation 
is extremely complex, and many efforts have been necessary to examine this problem in the light 
of theoretical physics tools, such as symmetries and their breaking [33j|34j. In this scenario a 
solution has been proposed by Parisi, with the well known Replica Symmetry Breaking scheme 
(RSB), both solving the SK-model by showing a peculiar "picture" of the organization of the 
underlaying microstructure of this complex system [28], as well as conferring a key role to the 
replica-trick method itself |39| . 

The physical relevance, and deep beauty, of the results, obtained in the frame of the replica 
trick, have prompted a wealth of further research, in particular toward the objective of develop- 
ing rigorous mathematical tool for the study of these problems. Let us recall, very schematically, 
some of the results obtained along these lines. Ergodic behavior has been confirmed in [l6j|23j, 
the lack of self-average for the order parameter has been shown in |35j, the existence of the 
thermodynamic limit in [22], the universality with respect to coupling's distribution in |15| . the 
correctness of the Parisi expression for the free energy in |21||38|. the critical behavior in [I], the 
constraints to the free overlap fluctuations in [2j[25j, and so much other contributions developed 
to give rise even to textbooks (see for instance |12j[18j[39j). 

Very recently, new investigations on ultrametricity started ([7J [8]) and allowed even strong state- 
ments dealing with the latter |30| . highlighting as a consequence the enquiry for techniques to 
prove the uniqueness of the Parisi solution, step by step. 

In this paper we match two other techniques, the broken replica symmetry bound |21| and the 
Hamilton- Jacobi method [20j[10j[19j, so to obtain a unified and stronger mathematical tool to 
work out free energies at various levels of RSB, whose properties are easily available as con- 
sequences of simple analogies with purely mechanical systems [19] , We stress that within this 
framework, the improvement of the free energy by increasing the replica symmetry breaking steps 
is transparent. 

In this first paper we show the method in full details, and pedagogically apply it for recovering 
the annealed and the replica symmetric solutions, then we work out the first level of RSB and 
show how to obtain the 1-RSB Parisi solution with its properties. 

The paper is organized as follows: In Section (2) the SK model is introduced together with its re- 
lated statistical mechanics definitions. In Section (3) the Broken Replica Mechanical Analogy is 
outlined in full details (minor calculations are reported in the Appendix), while Sections (4,5,6) 
are respectively dedicated to the annealed, the replica symmetric and the 1-RSB solutions of the 
SK model with our approach. Section (7) deals with the properties of the solutions and Section 
(8) is left for outlooks and conclusions. 
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2 The Sherrington-Kirkpatrick mean field spin glass 



The generic configuration of the Sherrington-Kirkpatrick model |26[ [27] is determined by the 
iV Ising variables a = ±1, i = 1,2, ... ,N. The Hamiltonian of the model, in some external 
magnetic field h, is 

H N (a,h;J) = — -= ^ Jij^aj - h ^ crj. (1) 

l<i<j<N l<i<N 

The first term in (pQ) is a long range random two body interaction, while the second represents 
the interaction of the spins with the magnetic field h. The external quenched disorder is given by 
the N(N — l)/2 independent and identically distributed random variables J^, defined for each 
pair of sites. For the sake of simplicity, denoting the average over this disorder by E, we assume 
each Jij to be a centered unit Gaussian with averages 

E(J fj ) = 0, E(4) = 1. 

For a given inverse temperatur43 /3, we introduce the disorder dependent partition function 
Zj^{(3,h; J), the quenched average of the free energy per site fpj{f5,h), the associated averaged 
normalized log-partition function aj^{(5,h), and the disorder dependent Boltzmann-Gibbs state 
to, according to the definitions 

Z N {/3,h;J) = ^2ex-p(-(3H N (a, h; J)), (2) 

<7 

-PMP,h) = jjElnZ N (p,h)=a N (P,h), (3) 
u(A) = Z N (/3, h; jy 1 '^2A(a)exp(-(3H N (a, h; J)), (4) 

where A is a generic function of a. 

Let us now introduce the important concept of replicas. Consider a generic number n of inde- 
pendent copies of the system, characterized by the spin configurations cr^\ . . . , a^ n \ distributed 
according to the product state 

n = JV xw (2) x ...xi>\ 

where each ui^ acts on the corresponding variables, and all are subject to the same sample J 
of the external disorder. These copies of the system are usually called real replicas, to distinguish 
them from those appearing in the replica trick |28| . which requires a limit towards zero number 
of replicas (n — > 0) at some stage. 

The overlap between two replicas a, b is defined according to 

^\^) = ^ E W> (5) 
Ki<N 



1 Here and in the following, we set the Boltzmann constant fca equal to one, so that /3 = 1/(A;bT) = 1/T. 
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and satisfies the obvious bounds 



1 < Qab < 1- 



For a generic smooth function A of the spin configurations on the n replicas, we define the 
averages (A) as 



where the Boltzmann-Gibbs average acts on the replicated a variables and E denotes, as usual, 
the average with respect to the quenched disorder J. 

3 Thermodynamics through a broken replica mechanical analogy 

Once introduced the model, let us briefly discuss the plan we are going to follow. 

In the broken replica symmetry bound (BRSB) [21 j it has been shown that the Parisi solution 

is a bound for the true free energy (the opposite bound has been achieved in |38|). This has 

been done by introducing a suitable recursive interpolating scheme that we are going to recall 

hereafter. 

In the Hamilton- Jacobi technique instead [20], it has been shown, by introducing a simple two 
parameter interpolating function, how to recover the replica symmetric solution trhough a me- 
chanical analogy, offering as a sideline a simple prescription, once the bridge to mechanics was 
achieved, to proof the uniqueness of the replica symmetric solution. 

The main result of this paper is that the two approaches can be merged such that even the 
recursive interpolating structure of the BRSB obeys a particular Hamilton- Jacobi description. 
This result has both theoretical and practical advantages: the former is a clear bridge among 
improving approximation of the free energy solution and increasing the levels of RSB, the latter 
is a completely autonomous mechanical tool by which obtain solutions at various RSB steps in 
further models. 

The task is however not trivial: the motion is no longer on a 1 + 1 Euclidean space-time as in 
|20| but lives in K + 1 dimensions such that momenta and mass matrix need to be introduced. 
To start showing the whole procedure, let us introduce the following Boltzmannfaktor 



where both the JijS and the Jfs are standard Gaussian random variables jV[0, 1] i.i.d. The t 
parameter and each of the x a maybe tuned in R + . We will use both the symbol x as well as 
(xi, ...,xk) to label the K interpolating real parameters coupling the one body interactions. K 
represents the dimensions, corresponding to the RSB steps in replica trick. Let us denote via E a 
each of the averages with respect to each of the J a 's and Eq the one with respect to the whole 
Jij random couplings. Through eq. (JJJ we are allowed to define the following partition function 




(6) 
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Zj<[{t; x\, . . . , xk) and, iteratively, all the other BRSB approximating functions for a = 0, . . . , K: 

Z K = Z N = Y,B({a}; X ,t), (8) 



Z™« x ee E a (Z™"), 



(9) 



Z™ 1 = E 1 {Z™ 1 ). 
We need further to introduce the following interpolating function 



and define, for a = 1, . . . , K , the random variables 

fa 

and the generalized states 



,xk) = j^E logZ , 



E a (Z^) 



(10) 



(11) 



(12) 



w o (0 = E a+ i ...E K (/ 0+ i . . . /j<-w(.)) , (13) 
the whole in complete analogy with the "broken prescriptions" [21] , 

Of course the corresponding replicated states £l a are immediately generalized with respect to 
each of the u a state introduced above. 

Overall, for a = 0, . . . , K, we further need the averages 

(14) 



0« = £(/l.../afi«(OJ- 

While it is clear that, when evaluated at t = f3 2 and x = 0, our interpolating function a(t, x) 
reproduces the definition of the quenched free energy, when evaluated at t = (which a proper 
choice for the x parameters that we are going to show), it reproduces the Parisi trial solution 
f(q = 0, y = h) at the given K level of RSB: 

a N (t = 0;xi, . . .,x K ) = 
1 



N 



log 



K 



m K -, 



E K exp J2^a~ 1 ^2 J i <Ji 



\a=l 



(15) 



Even though far from being trivial, this is an essential feature of mean field behavior even in the 
disordered framework; in fact, in the thermodynamic limit the connected correlation inside pure 
states should go to zero bridging the two body problem to a (collection of) one body model, or 
better "high temperature model", whose partition function factorizes: 



E 



exp ( v^Y. J >') = 2 " II cosh ( E ^ J n ■ 

\a=l i / i \a=l / 



(16) 
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such that, averaging over the , we get 



e k (y ex P {y^ ) = 

= 2 Nm « H I dfi(z K ) cosh" 1 * (E + 



ZK^JXK , 



(17) 



and so on. Even taking the external field h, which is again encoded in a single body interaction 
and is simply added into the hyperbolic cosine, we get 

a N (t = 0;xi, . . . ,x K ) = log 2 + 



+ loe 



j d^zx)... J dfi(z K ) cosh mA " (^^y/x~z a + j3h 



In the case where x a = (3 2 {q a — q a -i) the second term does coincide sharply with the solution of 
the Parisi equation |28| . 

Let us now define S(t, x) as Principal Hamilton Function (PHF) for our problem: 



S(t;xi, . . .,x K ) = 2 a(t;xi,. ..,x K 



a=l / 



(18) 



As proved in the Appendix, the (x, t)-streaming of S(t;xi, . . . ,xk) are then 

1 K 

d t S(t;xx,...,x K ) = -- ^2(m a+ i -m a )(ql 2 ) a , 

a=0 
1 K 

d a S(t;xi,...,x K ) = --^2(m b+1 -m b ){q 1 

b=a 



112 b- 



(19) 
(20) 



It is then possible to introduce an Hamilton- Jacobi structure for S(x,t), which implicity defines 
a potential V(t;x\, . . . ,xk), so to write 



1 K 

d t S(t, x) + - E d a S (M^) ab 8 b S + V(t, x) = 0. 



(21) 



a,b=l 
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The kinetic term reads off as 

K 



T = - d a S(t;x 1 , . . . ,x K ) (M 1 ) ab d b S(t; x x , . . . , x K ) 
o,6=l 

^ K K K 

= 2 { M ^)ab y^y^(m e+ i - m c ){q 12 }c(m d+1 - m d )(qi 2 ) d 

a,b=l c>a d>b 

1 K 

9 D cd(m c+ i - m c )(qi2)c(m d+ i -m d )(qi 2 }d, 



2 

c,d=l 



where we defined 



(22) 



c d 



a=l 6=1 

By the inversion of the mass matrix 



D cd (m c+ i - m c ) = 5 cd (24) 



we obtain the expression 
T 



~<T(mc+i-m c )(qi2} 2 c (25) 

c=l 

1 K 1 

= 2^( m c+i-^c){q\2)l- ^{rni-rnQ){qi 2 )l. (26) 
c=0 

Condition (|24p determines the elements of the inverse of the mass matrix M -1 . 
In particular we stress that it is symmetric and the not zero values are only on the diagonal and 
all of them respecting (M _1 ) aiC1+ i = (M _1 ) a+ i >a : 

(M" 1 )!! = , (27) 

7712 — mi 

(M _1 ) aa = + , (28) 

m a +i -m a m a - m a -i 

{M-^a+x = , (29) 

ma+i - m a 

all the others being zero. 

The elements of the mass matrix M are determined by the equation 

Y J M ab{M- x ) bc = 5 ac , (30) 
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and it is immediate to verify that the following representation holds: 

M ab = 1 - m (aA6 ) . (31) 

With this expression for the matrix elements, by substituting eq.s (|19p and (|26p into (j2Tj) we 
obtain the expression for the potential such that overall 

1 K 

dtS(t;x h ...,x K ) + -Y J d aS(M- 1 ) ab d b S + V(t;x 1 ,...,x K ) = 0, (32) 

a,b=l 

1 K 1 

V(t;x 1 ,. . . ,x K ) = ^^{m a+1 -m a )({ql 2 ) a - {q 12 )l) + -{m 1 -m ){qi 2 )l. 

a=0 

Once the mechanical analogy is built, it is however prohibitive solving the problem as it is (i.e. 
integrate the equations of motion) ; instead we propose an iterative scheme that mirrors the replica 
symmetry breaking one: at first, by choosing K = 1, we solve the free field solution (we impose 
V(t,x) = 0) and we recover the annealed expression for the free energy. This is consistent with 
neglecting the potential as it turns out to be the squared overlap. Then, we avoid perturbation 
scheme to deal with the source but we enlarge our Euclidean space by considering K = 2. 
Again we work out the free field solution to obtain the replica-symmetric expression for the free 
energy, consistently with neglecting the potential; in fact the source we avoid, this time, is the 
variance of the overlap: a much better approximation with respect to K = 1. We go further 
explicitly by considering the K = 3 case and we get the 1-RSB solution in the same way (and 
so on). Interesting we discover that there is a one to one connection among the steps of replica 
symmetry breaking in replica trick and the Euclidean dimension in the broken replica mechanical 
analogy. The latter however incorporates, in a single scheme, even the annealed and the replica 
symmetry solutions. 



4 K = 1, Annealed free energy 

Let us now recover some properties of disordered thermodynamics by studying the K = 1 case 
so to show how the solution of the free problem coincides with the annealed expression. 
We assume x{q) = m\ = 1 in the whole interval [0, 1]. 

We show now that, within our approach, this implies a reduction in the degrees of freedom where 
the Hamilton- Jacobi action lives, such that the PHF depends by t only. 
The dynamics involves a 1 + 1 Euclidean space-time such that 

Z 1 = Z K = Z N = Y^ exp ^ft/NH N (a; J) + V^Y, Ji(J ^j ' ( 33 ) 
Zq is consequently given by 

Z = E X Z X = exp (jxj ex P (Vt/NH N (a; J)) . (34) 
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This implies, into the interpolating function, a linear and separate dependence by the x 

a(t, x) = ^ + ^E log exp (^ft/NH N {a; J)) . (35) 



N 

(T 

The x-derivative of a(t, x) is immediate, while for the i-one we can use the general expression 
previously obtained (cfr. eq.s (|19l20p ) 

dta = ~[l-(«? a )o], (36) 
d x a = -. (37) 
As a straightforward but interesting consequence, PHF does not depend on x and we get 
S(t,x) = 2a(t,x)—x — - 

= ^E log exp [y/tjNH N (a- J)) - |, (38) 

dtS = ~(q 2 12 ) , (39) 

8 X S = v(t,x) = 0, (40) 

(41) 

where v(t) defines the velocity field, which is identically zero such that x(t) = xq. 

In this simplest case, the potential is trivially the i-derivative of S(t, x) with a change in the 

sign, (the averaged squared overlap): 

V(t,x) = ±{q> 2 ) . (42) 

Now we want to deal with the solution of the statistical mechanics problem. As we neglect 
the source (we are imposing (q , f 2 )o = 0), we can take the initial value for S(x,t) as it must be 
constant overall the space-time. 

S = 5(0) = 2 log 2, (43) 
and, consequently, we can write the solution of the problem as 

a{t,x) = log2 + | + ^. (44) 

At this point it is straightforward to obtain statistical mechanics by posing t = (3 2 and x = 0: 

a N {P) = log2 + ^, (45) 
which is exactly the annealed free energy. 
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5 K = 2, Replica symmetric free energy 



In this section, by adding another degree of freedom to our mechanical analogy, we want to 
reproduce the replica symmetric solution of the statistical mechanics problem. 
We deal with K = 2. The order parameter is now taken as 

So 



Qi = Q, q2 = qx = l 

mo = mi = 0, m 2 = mx = 1, rris = rriK+i = 1- 



(47) 
(48) 



The auxiliary partition function depends on t and on the two spatial coordinates x\ and x 2 : 



Z N (t;xi,x 2 ) = ^expi y/t/NH N (a; J) + y/xi^JiVi + y/x^^jfai J , 

a \ i i / 

and with the latter, recursively, we obtain Zq. 

Zr = Z 2 = Zn, 
Z 1 = (E 2 Z™ 2 )^ = E 2 Z 2 , 
Z = (EiZ™ 1 )^ . 

The function Z\ can be immediately evaluated by standard Gaussian integration as 



Z\ = exp 

Concerning the function Zq we can write 
{E x Z'^)^l = exp 
= exp 



— log Ei [exp (mi log Z x )\ 
m\ 

— log Ex \l + mi log Zi + o(m?)l 
mi 



1 



— [mitfi log Z x + o(m?)] 
mi L J 



= exp 
= exp Ei log Z x + o(mi), 



(49) 



(50) 
(51) 
(52) 



(53) 



and consequently 



Zq = exp Ei log Zi . 



(54) 
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In this case, our interpolating function reads off as 



X2 1 

a(t, xi,x 2 ) = — + jzEqEi log 



(55) 



Again by using the general formulas sketched in the first section (cfr. eq.s (|19l20p ) we get for 
the derivatives 



d t a 
d Xl a 



d X2 a 



[1 " <9i 2 2 >i] , 
[1 " (9i2)i] , 



Evaluating our function at t = 0,xi = x®, x 2 = x\ we easily find 



Let us introduce now the K = 2 PHF 

5(t;xi,x 2 ) = 

together with its derivatives 

d t S = 

d xl S = 
d X2 S = 



X\ Xo t 

2 [a--±- -f - - 



1 



(9i 2 ) 1. 



2 
0. 



-(912) 



(56) 
(57) 
(58) 

(59) 
(60) 



(61) 

(62) 
(63) 



We observe that, even in this case, there is no true dependence by one of the spatial variables 
(x 2 ): this is due to the constant value of the last interval txik = fni where the order parameter 
equals one and can be Gaussian-integrated out immediately into the corresponding Z 2 getting 
the pre- factor exp^iVxij). 

As a consequence, we can forget the mass matrix as there is no true multidimensional space. 
Let us write down the Hamilton- Jacobi equation 



d t S(t, Xl ) + ^ (d Xl S(t,xi)) 2 + V(t, Xl ) = 0. 



(64) 



The potential is given by the function 



V(t,x 1 ) = -(( q 2 12 ) 1 -( qi2 )l) 



(65) 
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where 

(g? 2 )i = E E 1 f 1 n 1 (q 2 12 ) = E E 1 f 1 ^ £ (E 2 f 2 u(a i a ] )) 2 . (66) 

ij 

When taking x\ = and t = f3 2 the variance of the overlap becomes the source of the streaming. 

n/5 2 ,0) = i((g? 2 )-(g 12 ) 2 ). (67) 
As usual in our framework, we kill the source (i.e. F(i,x) = 0), and obtain for the velocity 

g(rr?) = -«i(0,x?) = ^ dji(;z) tanh 2 (zy^?). (68) 

This is the well known self-consistency relation of Sherrington and Kirkpatrick, namely 

q(P) = J dn(z) tajah 2 (py/Qz). (69) 

The free field solution of the Hamilton-Jacobi equation is then the solution in a particular point 
(and of course the choice is 5(0, x®) which requires only a one-body evaluation) plus the integral of 
the Lagrangian over the time (which is trivially built by the kinetic term alone when considering 
free propagation). Overall the solution reads off as 

S(t,x 1 ) = S(0,x° 1 ) + ±?(x° 1 )t, (70) 
by which statistical mechanics is recovered as usual, obtaining for the pressure 

a{t;x u x 2 ) = log2 + J dfi(z) log cosh(^z) + |(1 - q) 2 + y , (71) 

that corresponds exactly to the replica-symmetric solution once evaluated at x\ = x 2 = and 
t = fi 2 and noticing that = x(t) = x\ — qt. 

Within our description it is not surprising that the replica symmetric solution is a better de- 
scription with respect to the annealing. In fact, while annealing is obtained neglecting the whole 
squared overlap (q 2 2 ) as a source term, the replica symmetric solution is obtained when neglect- 
ing only its variance. 

Of course, nor the former neither the latter may correspond to the true solution. However we are 
understanding that increasing the Euclidean dimensions (the RSB steps in replica framework) 
corresponds to lessening the potential in the Hamilton-Jacobi framework, and consequently re- 
ducing the error of the free field approximation toward the true solution. 

6 K = 3, 1-RSB free energy 

The simplest expression of x(q) which breaks replica symmetry is obtainable when considering 
K = 3, 

= q < qi < q 2 < q 3 = 1, (72) 
= mi < m 2 = m < 7713 = 1. (73) 
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With this choice for the parametrization of x(q) the solution of the Parisi equation 



V + ifyf + - 2 <dyf? = (74) 

is given by 

f(0,h;x,P) = ^ J d(i(zi) log J dfi{z 2 ) cos\\ m [f3(^zi + ^q 2 - q\z 2 + h)] + 

+ ^(l-ffi), (75) 

and, using a label P to empathize that we are considering Parisi prescription, the pressure 
becomes 



a P {fi,h;x) = log 2 + f(0,h;x^)-^ 2 j\x(q) 



dq 



+ 



log 2 - -/3 2 [(m - l)ql - 1 - mq\ + 2g 2 ] + 

— y log y d//(z 2 ) cosh m [/3(^/gl2:i + V<?2 - q\Z2 + h)]. 



(76) 



Now we want to see how it is possible to obtain this solution by analyzing the geodetics of 
our free mechanical propagation in 3 + 1 dimensions. 



Let us define 



Z N (t;x 1 ,x 2 ,x 3 ) = ^exp 



^t/NH N (a; J) + ^ y 7 ^ E J ? a * 



o=l 



by which 



Z 3 = Zk = Zn, 

{Nx 3 \ x 
Z 2 = E 3 Z 3 = exp I J 2^ ex P 

Z x = {E 2 Z™) l ' m , 

Z = {E 1 Z™ 1 ) Vmi = exp(E 1 log Z 1 ) = exp 
For the interpolating function we get in this way 



a N (t;xi,x 2 ,x 3 ) = —E logZ = 



^t/NH N (a; J) + £ Jx~ a £ Jfo 



a=l 



-El log£ 2 Z 2 m 
m 



(77) 

(78) 
(79) 

(80) 
(81) 

(82) 



f + ^o^log U 



exp y/i/NH N (a; J) + ^ v 7 ^ ^ 



a=l 
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while for the derivatives we can use the general formulas, so to obtain 



d t a 
did. 
d 2 a 
d 3 a 



-[l-m(^ 2 )i-(l-m)(^ 2 ) 2 ], 

^[1 - m(q 12 )i - (1 - m)(q 12 )2], 

i[l-(l-m)<gi 2 > 2 ], 
1 

2' 



Then we need to evaluate the interpolating function at the starting time: 



X3 



a N (0-x»,4,x») = ^+log2 + 



(83) 
(84) 
(85) 
(86) 

(87) 



+ — / d/x(zi)log 



dfj>(z 2 ) cosh r 



xj^i + \jx^z 2 



The K = 3 PHF, as usual and previously explained for the K = 1, 2 cases, does not depend 
on the last coordinate (i.e. X3), such that we can ignore it when studying the properties of the 
solution. 



S(t;x 1 ,x 2 ) 



-EqEi log E 2 



exp VtN/2K(a) + £ ^ E J ^ 



a=l 



iVm 

- xi - x 2 - i/2. 

and the derivatives, implicitly defining the momenta (labeled by Pi,p 2 ), are given by 



d t S 
d!S 



m 2 1 — m 
-y(?i2>i " -3- 

xi, x 2 ) = -m{q 12 )i - (1 - m){q 12 ) 2 , 



{112)2, 



d 2 S = p 2 (t;x!,x 2 ) = -(1 - m)(q 12 ) 2 . 
The kinetic energy consequently turns out to be 

1 — m, 

■(912)2. 



m , , 2 1 — m 2 



T = T^ + 2 

and the potential, which we are going to neglect as usual, is given by 

V(t; Xl ,x 2 ) = \[m «<&)! - < gi2 )?) + (1 - m) «g? 2 ) 2 - <gi 2 >|)] • 
By having two spatial degrees of freedom, the mass matrix has a 2 x 2 structure now 

u ( 1 "' 1 '" 

-1/m l/[m(l — m)\ 



(88) 



(89) 

(90) 
(91) 



(92) 



(93) 



(94) 
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M 



1 1 — m 
1 — m 1 — m 



(95) 



Note that the eigenvalues of the mass matrix are always positive defined for m £ [0, 1]. 
We can determine now the velocity field 



vi(t;xi,x 2 ) = ^2{M )ibPb = -(qi2)i, 
6=1 
2 

v 2 (t;x 1 ,x 2 ) = ^2(M~ 1 ) 2b p b = (q 12 ) 1 - (q 12 } 2 . 



(96) 
(97) 



6=1 



So we get all the ingredients for studying the free field solution (the one we get neglecting 
the source). In this case the equations of motion are 



zi(t) = zS-(gi2>i(0;s?,x!j)t = x?-git 

x 2 (t) = x^ + (q 1 -{q 12 )i(0;x^x^))t = x^ + (q 1 -q 2 )t 



(98) 
(99) 



and we can see that q\ and q 2 satisfy the self-consistency relations in agreement with the replica 
trick predictions 



qi 

D{z) 



j dn{z) 

j dn{z) 
P(y/qlz + ^q 2 - qiy), 
j d/i(y) cosh" 1 0(z,y). 



D- l {z) J dfi(y) cosh m 6 (z,y)tanh 6 (z,y) 
D- l {z) [ d^(y)cosh m e(z,y)tanh 2 e(z,y) 



(100) 

(101) 
(102) 
(103) 



The PHF is obtained in coherence with the previous cases and obeys 

S(t; Xl ,x 2 ) = S(0; s?, x%) + T(0; x?, x° 2 )t, 

by which 

X\ x 2 t 



a(t;x 1 ,x 2 ,x 3 ) - 
and, remembering that 



2 4 



a(0;x?,x0,^)-|-| + T(0;x?,x°)^ 



X\ — x\ 



-Qit, 

x 2 -x\ = (<?i - g 2 )i, 



(104) 



(105) 



(106) 
(107) 
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we get the thermodynamic pressure in the space-time coordinates: 



a(t;x 1 ,x 2 ,x 3 ) 



t 



9 + log 2 - ^[-1 + 2q 2 - mqf - (1 - m)q 2 2 ] 



(108) 



+ 



m 



d\i{z 2 ) cosh r ' 



x\z\ + \]x%z 2 ) 



In order to get the statistical mechanics result, as usual, we need to evaluate the latter in t = /3 2 , 
x i = x 2 = %3 = 0, from which x\ = q~\t and x 2 = (q~i — q2)t, gaining once again (|76p . 



7 Properties of the K = 1, 2, 3 free energies 

In the previous sections, we obtained solutions for the Hamilton- Jacobi equation in the K = 1, 2, 3 
cases, without saying anything about uniqueness. For K = 1, the annealed case, there is no true 
motion so it is clear that there is just a single straight trajectory, identified by the initial point 
xq = x, intersecting the generic point (x,t), with x,t > 0. 

In the K = 2 problem, well studied in |20j, one can show uniqueness by observing that 
the function t(xo), representing the point at which the trajectory intersects the x-axis, is a 
monotone increasing one of the initial point xq, so that given x, t > 0, there is a unique point xq 
(and velocity q(xQ), of course) from which the trajectory starts. 

For K = 3, the problem becomes much complicated, because we now have to consider motion 
in a three dimensional Euclidean space, proving that given the generic point (xi,x 2 ,t), with 
x\ > 0, x 2 > 0, t > 0, there exists a unique line passing in (xi,x 2 ) at time t. 

So let us consider the functions 



F(xi,t;x®,x 2 ) = x x -xi +q 1 (x1,x < ^)t, 



G{x 2 ,t 



= x 2 - x 2 + q 2 (xi,X2)t 



9i(x?,x§)t. 



(109) 
(110) 



These functions vanish in the points corresponding to the solutions of the equations of motion, 
and in particular for all the At = (xi = 0, x 2 = 0,t > 0; x\ = Q,x 2 = 0). Labeling with d\ 
and d 2 the partial derivatives with respect to x\ and x®, the Dini prescription tells us that if the 
determinant of the Hessian matrix 



d(F,G) 
d(x^,x^) 



diF d 2 F 
d x G d 2 G 



(111) 



is different from zero in a neighborhood of A tl then we can explicitate and x® as functions of 
Xi, x 2 and t, in such neighborhood. This means that the initial point and the velocities, which 
depend on it, are univocally determined by X\ , x 2 and t via the equations of motion. 
Calculating the determinant we find 



f, o o» = (- 1 + 5 i9i*)(- 1 + a 2g2t-9 2 git)-(a 2 g 1 t)(+9 1 g 2 t-9 1 g 1 t) 
o{x 1 ,x 2 ) 

= 1 + {d 2 qi - diqi - d 2 q 2 )t + {diqid 2 q 2 - d 2 qidiq 2 )t 2 . 



;ii2) 
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so we should ask, for all > and x\ > 

A = (d 2 qi - diqi - d 2 q 2 ) 2 - 4(<9i<7i<9 2 cz 2 ~ d 2 qidiq 2 ) (113) 
to be negative, or in case A > 0, the zeros of 



t± _ ~d 2 qi + glgl + d 2 q2 ± VA 
2(d 1 q 1 d 2 q 2 - d 2 qidiq 2 ) 

correspond to non-invertibility points. 

The expression we obtain for the determinant is quite untractable, however we can show 
uniqueness in a neighborhood of the initial point x\ = 0, x^ = 0. The motion starting from this 
point has zero velocity, and we saw that it gives the high temperature solution for the mean field 
spin glass model: Remembering that the transition to the low temperature is continuous, we can 
expand the Hessian for small values of x\ and x® and observe that, for x\ = 0, x 2 = 0, t = f3 2 , 
the equations of motions become 

x? = fqi (115) 

A = /3 2 fe-gi). (116) 

When x\ — >• and x® — > we have also q\ — >• and q 2 — > 0, so we have an expansion close to the 
critical point (which is the only region where the control of the unstable 1-RSB solution makes 
sense for the SK, being the latter oo-RSB). 

For qi and q 2 we have, retaining terms until the second order: 

?i(xS,x§) ps x? - 2(1 - m)x° 1 x° 2 - 2(x\f (117) 
f x\ + x° 2 + mx 2 {xl + 2x1) (118) 



q 2 (x°,x 2 ) 



and consequently 



9iqi(x?,x^) ps 1 - 2(1 - m)x\ - Ax\ (119) 

3 2 q x (x\,x%) ps -2(1 -m)x? (120) 

01<fc(ar?,a:§) ps 1 + 2/71x1] (121) 

a 2 ^2(x?,x^ ps l + 2mx? + 2mx^. (122) 



Substituting in fjl 121) we find 

|^fy - l-2[l-x?- ( l-2m)x°]t 

+ [l - (2 - m)x\ - 2(1 - 2m)x% + 2m(m - 4)x?x 2 + 

-%m{x\f -Ara{l-m)(xl) 2 ]t 2 . (123) 

and, for x^x^ = (which corresponds to expand the velocities up to the first order in x^ and 
x®) we simply obtain 

a(xp Xg) 

This means that in a neighborhood of x\ = x 2 = we have uniqueness, provided that we are 
not exactly at the critical point t = j3 2 = 1. 
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8 Outlooks and conclusions 

In this paper we enlarged the previously investigated Hamilton- Jacobi structure for free energy 
in thermodynamics of complex systems (tested on the paradigmatic SK model) by merging this 
approach with the Broken Replica Symmetry Bound technique. At the mathematical level the 
main achievement is the development of a new method which is autonomously able to give the 
various steps of replica symmetry breaking (of the replica trick counterpart). At a physical 
level this methods clearly highlights why increasing the steps of RSB improves the obtained 
thermodynamics mirroring these increments in diminishing the approximation of a free field 
propagation in an Euclidean space time of an enlarged free energy, which recovers the proper 
one of statistical mechanics as a particular, well defined, limit: the main achievement is paving 
an alternative way to understand RSB phenomenon. 

However, when increasing the steps of RSB (making smaller the potential we neglect, and so 
smaller the error) there is a price to pay: each step of replica symmetry breaking enlarges by 
one dimension the space for the motion of the mechanical action. As a consequence the full RSB 
theory should live on an Hilbert space: this still deserve more analysis, however the method 
is already clear and several application may now stem: for example P-spin above the Gardner 
critical temperature could be solved exactly as well as a consistent part of the plethora of models 
born in disordered system statistical mechanics whose solutions implies only one step of RSB. 
We deserve to investigate both the K — > oo limit to complete the theory as well as its simpler 
immediate applications. 

Appendix: Streaming of the interpolating function <S(£,x) 

In this section we show in all details how to get the streaming of the interpolating function (|lip 
The i-streaming of the interpolating function a(t, x) is given by the following formula: 




(125) 



To get this result, let us start by 




(126) 



and, as it is straightforward to show that 




(127) 



by iteration we get 



Z \ X , t)d t Z ( X , t) = E 1 ... E K (h . . . f K Zx\ X , t)dtZ K (x, t)). 



(128) 
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The i-derivative of Zk is then given by 

Z- X (x, t)d t Z K (x, t) = — = JiM^j), (129) 

ij 

from which 

Z \ X , t)d t Z (x, t) = —L= f K JiM°i°i)) , (130) 

where we labeled with i£ the global average overall the random variables as there is no danger 
of confusion. All the terms into the sum can be worked out integrating by parts: 

K 

E(h... f K Jij^{(yi<yj)) = ^2 E (/i . . . dj..f a . . . f K u{aiaj)) 

a=l 

+ E(f 1 ...f K dj iJ u>(* i * j )). (131) 

So we need to calculate the explicit expression of the derivatives with respect to Jij of both f a 
as well as u(aiGj). For the latter, it is easy to check that 

dj^icxicxj) = ^ (1 - o/Wi)) , (132) 

while for the / a 's we have 

d. h] f a = m a f a (Z-^x, t)d Ji3 Z a (^ t)) - m a f a E a f a {Z~\^ t)dj tj Z a (x, t)) . (133) 
By using the analogy of (|127[) we get 



Z- 1 (x, t)a Jy Z (x, i) = K+i • • • Ex (f a+1 . . . f K z~ K x d Ji3 Z K ) 

N 



UJaiViCTj), (134) 



such that 



dj l3 fa = m a f a J — (oj a (ai<7j) - w a -i(<7i<7,-)) • (135) 



Substituting (|132p and (|135p into f j 1 3 1 P we obtain 

K 



E(fl... f K JijU{(Ti<Jj)) = J— [E (/l . . . faOJaiViCTj) . . . f K u(cTiCrj)) 

a=l 

- E(fx-.. f a -lU} a -l((Ji(Jj) ■ ■ ■ fKu(<Ji(Tj))) 

+ ^E(f 1 ...f K (l-u\a i a j ))). 



(136) 
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Overall, an explicit expression for the eq. f|126[) is given by 
1 K 

9tSl = 4/V2^^ ma ^° ' ' ' Ea ^ ' ' ' faUa(cnvj)E a+1 . . . E K f a+ i . . . f K uj(aiaj) 

a=l ij 

- E ... £ a _i/i . . . f a _ x u a -\(?i<Jj)E a . . . E K f a . . . f K u)(aiaj)] 
^.../^(l-w 2 ^)). 



4AT2 



(137) 



ij 



Once introduced the overlap, we can write the result: 
K 



a=l 

1 ( K K \ 

= 4 S ma ^l2>a -^2m a+1 {qj 2 ) a + m K+1 (ql 2 ) K + 1 - {q\ 2 ) K 

\a=l a=0 / 



(135 



Now let us focus on the x-streaming of the interpolating function a(t,x) and show that it is 
given by the following formula: 



d a a N (x, i) = ~ ^1 - y2(m b+1 - m b )(q 12 (x, t))^ 



In analogy with the i-streaming we have 



d a a N (-K,t) = —E Z 1 (-xL,t)d a Z (-x,t), 

Z 6 _1 (x,t)5 a ^ 6 (x,t) = E b+1 (f b+1 Z b ^ 1 (x,t)d a Z b+1 (x,t)) , 

>Z- 1 (x,t)a a Z (x,t) = E 1 ...E K (f 1 ...f K Z-\x,t)d a Z K (x,t)), 

Zx\x,t)d a Z K (yL,t) = — ^^JfcD^), 



by which 



d a a 



1 1 



iV2^ 

Again by integrating by parts we have 



^2Eif 1 ...f K J?£,(a i )). 



(139) 

(140) 

(141) 
(142) 

(143) 
(144) 



1 1 N 
d a a = — tt^= 



N2J¥~ a 



i=l 



K 



Y,E{fi... dj«f b . . . f K £o{c?i)) +E (/!... /A-SjxiwCffi))] . (145) 



6=1 
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Let us work out the Jf by remembering that Zb's, and consequently f b s, do not depend on 

T b+1 T K 
i J • • • i J i ■ 



dj?h = { 



( 

m, 



if a > b 

■ a f a (V(x, t)dj«Z a (x, t)) if a = b (146) 

k m b f b (Z b ~ 1 (x, t)dj*Z b (x, t)) - m h f b E h f b (Z b -1 (x, t)dj»Z b (x, t)) if a < b. 



The same recursion relationship holds in this case as well: 

Z-\yL, t)dj«Z b (x, t) = E b+ \ ...E K (f b+1 . . . f K Z K \^ t)dj«Z K (x, t)) 

Furthermore 



Z K 1 (x,t)d J ^Z K (x,t) = y/x^uj{Gi 



from which we get 



Z b 1 (x,t)0jaZ b (x,t) = y/x^U} b ((Ti). 

Consequently, eq.s (|146p can be written as 

[ if a > b 

dj?fb = I m a f ay /x^L) a (<Ji) if a = b 

{ y/x^m b f b (w b ((Ti) - d) 6 _i(crj)) if a < b. 

The last thing missing is evaluating the derivative of the state 

djo-u}(ai) = y/x^ (l - u?{(Ji)) , 



(147) 

(148) 
(149) 

(150) 
(151) 



so to write, via the overlap, the analogous for the generalized states. Substituting eq.s f 1 1 5 1 () and 
(|150p . once expressed via overlaps, into (|145p we obtain eq. p39D . 
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